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Two-Kaon Correlations in Central Pb + Pb Collisions at 158 A GeV/c 
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Two-particle interferometry of positive kaons is studied in Pb + Pb collisions at 
mean transverse momenta < pt >~ 0.25 and 0.91 GeV/c. A three-dimensional 
analysis was applied to the lower px data, while a two-dimensional analysis was 
used for the higher px data. We find that the source-size parameters are consistent 
with the niT scaling curve observed in pion-correlation measurements in the same 
collisions, and that the duration time of kaon emission is consistent with zero within 
the experimental sensitivity. 

PACS Numbers: 25.75.Gz 



Experimental studies of high-energy nuclear collisions 
at the BNL-AGS and CERN-SPS accelerators (beam en- 
ergies from 10 to 200 GeV/nucleon) have revealed in- 
teresting features of hot and dense nuclear matter, and 
some characteristic signatures of a quark-gluon-plasma 
(QGP) phase have been reported |jl|. If the hadronic 
source is formed in a first-order phase transition from a 
QGP phase in the course of collision, the hadronic expan- 
sion may slow down due to a softening of the equation 
of state. In such a case, a long duration time of par- 
ticle emission is anticipated. Since a finite duration of 
particle emission increases the effective source size in the 
direction of particle velocity, and since the shape of the 
two-particle correlation function is related to the effec- 
tive source size, a difference between the widths of the 
peaks in the correlation functions in the direction of pair 
("outward") and perpendicular to it ("sideward") might 
be a signature of QGP g]. At SPS energies, systematic 
studies of particle correlations have been performed from 
p + Ato Pb + Pb collisions y-QI • From the pion correla- 
tion studies in the Pb -I- Pb collisions, NA44 [|| reported 
that the two transverse radius parameters in the out- 
ward and sideward directions in the longitudinal center 



of mass system (LCMS) frame (see below) are similar, 
implying no long duration time of emission. WA98 0] 
measured the correlation function with the generalized 
Yano-Koonin parametrization and found the Rq parame- 
ter, which reflects the duration of emission, is compatible 
with zero. While NA49 H observed a finite _Ro parame- 
ter in the Yano-Koonin-Podgoretskii parametrization to 
be approximately 3-4 fm/c, the value is small and not 
consistent with what would be expected from a strong 
first-order phase transition. All the experimental data 
support no long-lived intermediate hadron-parton mixed 
phase during the pion emission. There are discussions, 
however, that the pion correlation functions might be dis- 
torted due to a large amount of decays from long-lived 
resonances, while kaon measurements can serve as a more 
sensitive probe of the space-time evolution |p|-[10| . The 
kaon duration time was measured in S -I- Pb collisions, 
but a long-lived mixed phase was also excluded by the 
data Q. In this letter, we present the first results of 
K^K^ correlations in central Pb -f Pb collisions at the 
SPS energy and extract the duration time of the kaon 
emission. 

The data were taken at the CERN-SPS with a 158 



GeV/c per nucleon lead- ion beam incident on a lead tar- 
get, using the NA44 spectrometer at two laboratory an- 
gles: 44 and 131 mrad with respect to the beam axis. 
A description of the NA44 apparatus can be found else- 
where nM- The lead beam was transported in vacuum 
up to a beam counter. The signal amplitude from the 
beam counter ensures a single lead ion impinging on the 
lead target at the entrance edge of the first dipole mag- 
net. Behind the target, there were two multiplicity de- 
tectors characterizing each event. A set of scintillator 
bars was used primarily to generate a centrality trigger 
by measuring charged secondaries. A silicon pad detector 
with a full coverage in the azimuthal direction and seg- 
mented into 512 pads measured the dN/drj distribution 
of charged particles in the pseudorapidity range between 
1.5 and 3.3. 

The secondary particles were transported to the track- 
ing section through another dipole magnet and three 
quadrupole magnets. The quadrupole magnets control 
the three-dimensional momentum acceptance of the spec- 
trometer. Two sets of quadrupole settings, referred to as 
the "horizontal" and "vertical" focus settings, were used. 
The horizontal (vertical) setting favors a wide momentum 
acceptance in the horizontal (vertical) space p^ (py), and 
reduces the acceptance Py (px)- In the discussion which 
follows, the direction with wide acceptance for a given 
mode (e.g., Px for the horizontal mode) is referred to as 
the "favored" direction and the direction with narrow 
acceptance (e.g., py for the horizontal mode) is called 
"unfavored". For the small angle measurements, colli- 
mators between the first dipole and the first quadrupole 
magnets further reduced the acceptance in the unfavored 
direction. These collimators reduced the number of par- 
ticles in the acceptance without reducing the acceptance 
in the favored direction. The momentum settings of the 
spectrometer were 6 GeV/c for the small angle (low px) 
measurements and 7.5 GeV/c for the large angle (high 
Pt) measurements. 

After the magnets, there were three tracking cham- 
bers (PC, SCI, SC2), three threshold-type gas Cherenkov 
counters (CI, C2, TIC) and three scintillator hodoscopes 
(H2, H3, H4). The pad chamber (PC) and the strip 
chambers (SC) provided precise hit position for the track- 
ing algorithm. Although their position information was 
less precise, the hodoscopes were also used in track find- 
ing. CI was the primary device for pion/kaon separation 
in the small angle measurements, while C2 and TIC were 
used for pion/kaon separation in the large angle measure- 
ments. Tracks were reconstructed by requiring hits in 
the tracking chambers and hodoscopes on straight lines. 
The three-dimensional momenta were calculated from 
the track information through a matrix, of which the 
elements were determined by a Monte Carlo simulation 
combined with the TURTLE code for the particle trans- 
portation through the magnetic field. Particle identifi- 
cation used time-of-flight information between the beam 



counter and the hodoscopes along with pulse height in- 
formation from the Cherenkov counters. Figure n shows 
an example of such a plot in the small angle measure- 
ments and also the particle-identification (PID) selection 
adopted in this analysis. Contaminations of pion and 
proton tracks in the kaon-identified pair-samples were 
evaluated to be 0.15% and 3.9%, respectively. Therefore, 
contaminations of tttt and pp pairs in the KK samples 
were negligibly small. 
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FIG. 1. The particle identification using the mass-squared 
and the CI pulse height. Kaon and proton peaks are seen 
near the CI pedestal (70 ch.), while a bump at around 600 
ch. on CI and Mass^ « is of pion events accumulated with 
a hardware scaled-down trigger. On the top plane, a scatter 
plot is shown with the PID selection adopted in this analysis. 

At the small angle, the spectrometer covers the rapid- 
ity range of 2.9-3.3 with a pT window below 0.6 GeV/c 
(< Pt >~ 0.25 GeV/c). The rapidity coverage in the 
large angle measurements is 2.4-2.9 with the pt window 
of 0.7-1.4 GeV/c (< pr >w 0.91 GeV/c). Note that the 
colliding system is symmetric and the central rapidity is 
at 2.9, so the rapidity coverage in the high and low pT 
measurements overlap well. The trigger required a lead 
ion in the beam counter, a large number of secondary 
particles in the scintillator bars, and at least two hits in 
hodoscopes H2 and H3. Off-line we required at least two 
kaons in each event and found the centrality of the lower 
Pt data set to be around 10% of the most central events 
while that of the high pT sample was 18%. 

The two-particle correlation function is defined by 



Crawipi,P2) = 



P2{pi,P2) _ Real{pi,p2) 
Pi{pi)Pi{P2) Back{pi,p2) 



(1) 



where the numerator is the joint probability of detect- 
ing two particles with momenta pi and p2, while the 
denominator is the product of the probabilities of de- 
tecting single particles with momenta pi and p2- The de- 
nominator was obtained by mixing two tracks picked up 
from two randomly selected different events. Ten mixed 
background pairs were generated for each real pair to re- 
duce statistical uncertainties. There still remain several 



effects, which affect the true correlation function. The 
limited acceptance of the spectrometer distorts the real 
two-particle spectrum, and the finite momentum resolu- 
tion smears the peak in the correlation function. The 
repulsive Coulomb force separates particles of a pair go- 
ing out at pi ~ p2 and strongly suppresses the yield at 
q = Pi — P2 ^ 0. The background two-particle spec- 
trum generated by the event mixing method still con- 
tains effects of two-particle correlations, since a particle 
in a real event is always accompanied by another particle 
nearby in phase space. Since the degree of these effects 
depends on the strength of particle correlations of inter- 
est in a source, these effects were corrected through a 
Monte Carlo (MC) based iteration method. A C2 func- 
tion (described in the next paragraph) was assumed for a 
given set of source-size parameters (i?'s and A), and MC 
events were generated taking the spectrometer response 
into account. The wave integration method |l^ was em- 
ployed to simulate Coulomb effects in the finite source 
volume. The MC events were analyzed using exactly the 
same procedure as was applied to the experimental data. 
The correction factors to the Craw function were evalu- 
ated by comparing the resultant MC correlation function 
with the given input C2 function. The source size pa- 
rameters were deduced from the experimental correlation 
function after applying the correction factors to the Craw 
function. The iteration ends when the extracted source- 
size parameters agree with the given parameters within 
an acceptable accuracy. The detailed procedure is de- 
scribed in the first article in Ref. 0] and more discussion 
can be found in Ref. |jl3|. 

After the PID selection and quality cuts, we have 
around 20 x 10'^ pairs in both the horizontal and vertical 
modes at the lower px, while 17 x 10'^ pairs in the hori- 
zontal mode at the higher px- A three-dimensional fit is 
applied to the lower pT data. The observed momenta of 
each kaon pair are transformed to the LCMS, in which 
the momentum sum of two particles along the beam axis 
is zero, Pzi +Pz2 = 0. The momentum variables of inter- 
est are defined in this system. The average momentum 
of the pair in the frame is kx = {pti + PT2)/2. Ql is 
a projection of the momentum difference onto the beam 
axis, and Qt is the momentum difference perpendicular 
to the beam axis. Qt is further divided into two compo- 
nents. Qto is the component of Qt along fc-r, and Qts 
is the component of Qt perpendicular to both fc^ and 
the beam axis. The momentum resolution of Ql, Qto, 
and Qts, including effects of multiple scattering in the 
target, were around 10, 20, and 20 MeV/c, respectively. 
The correlation functions C2 in both the horizontal and 
vertical modes are simultaneously fitted using the maxi- 
mum likelihood method with 



C2{QTO,QTS,QL)^D{l + Xe- 



0, 



(2) 



where A is a factor introduced to express chaoticity 
of quantum states of the source and i?'s are vari- 
ables representing multidimensional radii of the system 
in question. Z? is a free parameter for normalization 
in each mode. A two-dimensional equation replacing 

—Q'^qR^q — Qts^ts ^^ ^1- H ^i*^ —QtRt is employed 
to fit the correlation function in the horizontal mode at 
the higher pT- The fit results are given in Table I, where 
the systematic uncertainties reflect the effect of (i) cut 
parameters to define a track, (ii) cut parameters to select 
pairs, (iii) momentum resolution, (iv) two-track resolu- 
tion, (v) momentum distribution of particle production 
in MC, and (vi) fitting to finite bins. Figure g shows 
projections of the correlation function onto each axis of 
the momentum difference where the projection is over 
the lowest 40 MeV/c of the other directions in the mo- 
mentum difference. In these plots, the solid lines show 
the results of fit projected in the same way as the data. 

TABLE I. Fit results of Gaussian parametrizations 
of the K^K^ correlation functions at low and high 
Pt- The errors are statistical and systematical ones. 



<PT> [GeV/c] 

A 

Rl [fm] 

Rt [fm] 

Rts [fm] 

Rto [fm] 

X^ld.o.J. 



0.25 



0.91 



0.84 ± 0.06 ± 0.07 0.61 ± 0.20 ± 0.16 
4.36 ± 0.33 ± 0.32 3.20 ± 0.54 ± 0.45 
n/a 3.59 ± 0.67 ± 0.97 

4.04 ± 0.28 ± 0.32 n/a 

4.12 ± 0.26 ± 0.31 n/a 

5139 /2978 117 / 107 
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FIG. 2. Projections of the K 
PT ~ 0.25 GeV/c (top) and pr 



K correlation functions at 
: 0.91 GeV/c (bottom). The 

solid lines show the projections of the fit with the Gaussian 

parametrization. 



The three source-size parameters {Rl, Rto, Rts) in 
the three-dimensional fit are quite similar to each other 



— as was observed in pion correlation measurements in 
the same colliding system [||. Since Rts and Rl repre- 
sent the geometric information of the source most di- 
rectly, they are compared in Fig. @ with those of pi- 
ous. The present data at itit ~ 0.55 GeV/c^ seems to 
stay on the niT scaling curve (dashed) which came from 
the pion correlation measurements. To confirm this ten- 
dency, we put Rt and R^ from the two-dimensional fit 
at the higher px {rriT ~ 1.0 GeV/c^) on the same plot. 
They also seem consistent with the scaling curve. A fit 
to the four data points in each plot with a single scaling 
curve, R = A/y^, gives A = 3.0 ± 0.2 fmGeV^/^c~i 
in both cases as shown with solid curves. 
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FIG. 3. The longitudinal and transversal source-size pa- 
rameters of K+K+ at pt ^ 0.25 GeV/c and 0.91 GeV/c 
(solid circle), compared with those of pions (open circle). The 
dashed curve shows a fit to _R = A/^/mr for pions, while the 
solid curve is the fit to the pion and kaon data-points in each 
plot. 

The dependences on rriT are predicted with models in- 
cluding hydrodynamical expansion in a source [HSI-Em. 
The experimental radius parameters are interpreted as a 
length of homogeneity which is in turn dependent on a ge- 
ometrical source size Rgeom and a thermal length Rtherm- 
The "boost invariant expansion" along the longitudinal 
direction leads to an expression Rl = T^Tf^lvfiTi ignor- 
ing the contribution from Rgeom- Assuming a freeze-out 
temperature Tq of 100-140 MeV, we could extract the 
freeze-out time r from the present A parameter to be 7~ 
10 fm/c, which is in good agreement with the WA98 ra] 
and NA49 M results. In a hydrodynamic model under 
certain conditions |l3] , one can derive an analytic expres- 
sion for the radii Rt — Rl — t^To/ijit- Our data are 
consistent with such a scenario. The common tot scaling 
for pions and kaons may imply that thermal freeze-out 
occurs simultaneously for both pion and kaons and that 
therefore they receive a common Lorentz boost. This is 
consistent with the hydrodynamic hypothesis. A similar 
conclusion can be drawn from the linear increase of the 
single particle inverse slopes with mass |16| . 

We derived the duration time At of kaon emission 
from quadratic difference of Rto and Rts in Eq. (H) 



as At — \/Rto ^ RtsIP-i where /3 is the transverse ve- 
locity of the kaon pair in the LCMS frame. We find 
At = 2.2 ± 5.2{stat) ± 6A{syst) fm/c. The kaon dura- 
tion time is short and similar to those observed for pions 
in the same colliding system and for kaons in the S 4- Pb 
collisions. The present result excludes simple scenarios 
of a prolonged mixed phase anticipated in a first-order 
phase transition from a QGP phase. 
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